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ABSTRACT.  Completing  preliniiniiry  work  initiated  by  \OTS  several  years  ago, 
photoelastic  tests  were  conducted  to  determine  the  maximum  tangential  stress  for 
a  family  of  star-grain  configurations  in  a  state  of  plane  stress.  Curves  are  pre¬ 
sented  that  give  the  value  of  the  ratio  r»f  tangential  stress  to  external  pressure  for 
web  fractions  from  20  to  00*^,  for  n/p  from  approximately  2  to  W'fw'  is  the  web 
thickness  of  the  grain,  and  p  is  the  ra<liiis  of  ilic  star  j)oint),  and  for  ,Y  (number  of 
star  points^  from  (1  to  8.  'I'he  accuracy  of  the  curves  is  estimated  to  be  within 
+  4®c.  The  effect  of  the  inverse  star  point  was  investigated  and  found  to  be  negli¬ 
gible  over  a  range  of  shapes  and  values  of  tin;  [>arameler  a'/a  (ratio  of  distances 
from  center  to  the  inverse  star  point  an<l  to  the  outer  star  point)  from  approximately 
zero  to  0.6,  The  experimental  residls  compared  favonibly  with  one  theoretical  cal¬ 
culation.  equations  are  given  for  extension  of  the  data  to  case-bondetl,  internally 
pressurized  grains  and  for  application  to  problems  of  tangential  strain.  An  appendix 
presents  an  empirical  equation  representing  the  data. 
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FOREWORD 

Some  preliminary  design  data  for  stresses  in  star  grains  were  developed  several  years  ago  by 
the  use  of  photoelastic  techniques.  Because  of  the  importance  stress  analysis  has  assumed  in  the 
development  of  current  and  proposed  propulsion  systems,  and  because  the  purely  theoretical  ap¬ 
proach  has  not  as  yet  developed  into  a  useful  engineering  tool,  extension  of  the  earlier  work  was 
deemed  desirable.  This  report  presents  the  results  of  the  additional  work,  which  increased  the 
number  of  parameters  and  improved  the  accuracy  of  the  results. 

The  work  was  performed  for  the  U.  S.  Naval  Ordnance  Test  Station  by  the  Mathematical 
Sciences  Corporation  (formerly  Southwest  Engineers)  under  Contract  No.  N123{60530)25914A 
during  the  period  November  1960  to  May  1961. 
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NOMENCLATURE 

a  Radius  to  outer  star  point,  in, 

_ u  Rudias  to  inverse  star  point,  in. _ 

-  T  -4.  lnver8e:atar:point^ 

^ --Radi  u^of  gt^inj^in - = 

rf  Differential 
D  Diameter  o^grwiaj  in, 

E  Yo^g* B  modulus 
/  Stress-optic  constant 
E  Fringe  constant  based  on  h  a  0.250  in. 

G  Shear  modulus 

Ki  Stress-concentration  factor  at  the  inner  surface  of  the  grain 
L  Load,  lb 

~n  Isochromntic  fringe  order 
N  Number  of  star  points 
Pg  Gage  pressure 
Pi  Internal  pressure 
Po  Kxternal  pressure 
p  Pressure  at  interface  of  grain  and  case 
P  Ratio  of  interface  pressure  to  internal  pressure 
w  Web  thickness  of  grain,  in. 


y  Shear  strain 
(  Strain 

i'  Poisson’s  ratio 

p  Radius  of  star  point 

a  Stress 
r  Shear  stress 


SUBSCRIPTS 

c  Pertaining  to  grain  case 
g  Gage 
i  Internal 
r  Radial 

z  In  a  direction  normal  to  the  plane  of  the  slice 
9  Circumferential 


iv 
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- ^INTRODUCTION 

The  extensive  use  of  large  solid-propellant  rockets  in  recent  years  has  intensified  the  need 
for  more  precise  methods  of  stress  analysis  of  such  systems.  The  complete  stress  analysis^  which 
is  but  a  part  of  the  real  task  of  determining  the  s^uctural  integrity  of  the  systemi  is  in  itself  an 
exceedingly  difficult  problem.  This  report  is  concerned  with  a  very  limited  aspect  of  the  general 
problem  and  deals  specifically  with  the  stress  and  strain  concentrations  that  exist  at  the  star 
points  of  a  long,  tubular  solid-propellant  grain  subjected  to  internal  pressure.  The  report  amplifies 
certain  preliminary  data  obtained  by  Southwest  Engineers  for  the  U.  S.  Naval  Ordnance  Test  Sta¬ 
tion  in  an  earlier  study  (Ref,  1),  which  provided  the  first  parametric  investigation  of  structural 
implications  induced  by  Imllistic  changes  in  the  charge.  Because  of  the  unusual  and  irregular 
shapes  of  rocket-grain  cross  sections,  it  is  usually  impractical  to  employ  analytical  stress- 
analysis  methods,  and  hence  alternate  techniques  become  attractive. 

Following  the  previous  investigation,  n  photoelastic  approach  was  used  to  improve  the  ac¬ 
curacy  of  the  preliininury  data  and  oxtend  the  applications  pointed  out  by  Ordahl  and  Williams 
(lief,  1).  The  tests  to  bo  discussed  concentrate  upon  a  family  of  star  grains  and  present  stress- 
concentration  data  believed  to  bo  accurate  within  ±i%  over  the  range  of  design  interest. 


STATEMENT  OF  PROBLEM  TO  BE  CONSIDERED 


1  ho  structural  integrity  of  solid-propelinnt  grains  is  a  problem  that  has  received  much  atten¬ 
tion  in  the  past  several  years  (Ref,  2),  I'lie  major  difficulties  are  due  to  the  need  for  a  viscoelastic 
analysis  and  the  lack  of  an  adequate  failure  criterion.  Advances  in  both  arcus  arc  hampered  by  the 
complicated  grain  geometry.  It  is  the  aim  of  this  report  to  present  one  set  of  data  showing  the  ef¬ 
fect  of  the  geometry  on  the  magnitude  of  the  maximum  stress  caused  by  a  pressure  loading.  While 
the  analysis  applies  specificiilly  to  an  elastic  body,  it  may  be  shown  (Ref.  2)  that  for  cases  where 
the  boundary  conditions  are  imposed  essentially  in  terms  of  stresses,  the  viscoelastic  and  elastic 
stress  distributions  are  identical.  Hence,  for  engineering  puqioses  it  is  found  llml  the  viscoelastic 
stress  distribution  in  the  cross  section  of  a  long  tubular  grain  may  be  obtained  by  photoelastic 
studies, 

Photoelastic  experiments  are  conducted  on  thin  slices  for  which  plane-stress  assumptions 
are  used,  i.e.,  the  stress  normal  to  the  slice,  (Tg,  is  taken  to  be  zero.  These  data  are  analytically 
transferable  to  a  long-tube  grain  where  the  axial  strain,  r*,  is  either  zero  or  constant.  This  is 
therefore  the  condition  to  which  the  data  will  apply;  specifically,  they  are  not  expected  to  apply 
near  the  ends  of  the  grai.*.  or  near  abrupt  changes  in  cross  section  such  as  at  a  head-  or  after-end 
closure. 

It  should  also  be  stated  that  the  rocket  grain  is  assumed  to  be  composed  of  an  isotropic,  ho¬ 
mogeneous  continuum,  notwithstanding  the  presence  of  filler  particles  and  possible  voids  in  the 
viscoelastic  binder  of  the  propellant. 

The.  geometry  of  the  specific  family  of  cross  sections  to  be  tested  is  shown  in  Fig,  1.  The 
maximum  stress,  as  shown  in  Fig.  2  by  the  concentration  of  the  photoelastic  fringes  around  each 
star  point,  occurs  at  the  star  point,  which  lies  along  a  radial  line  where  the  shear  stress,  is 
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f^t  rodiiis  of  grnint  <t,  rodius  to  outor  «tar  point;  a\  riidioi?  to  invorse  star 
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FIG.  2.  Five-Point-Star  Grain  Subjected  to  External  Pressure. 
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zero  by  syrainetry.  Hence,  the  radial  and-circHmferRiitial-s&osjies,  (7^  and  a respectively^are 
the  principal  stresses  that  permit  immediate  application  of  the  fundamental  stress-optic  law 


Oi-02^^fn/h  (1) 

where 

Olj2  “  the  principal  sfiessea 

/=  8tress»optic  con8tant  for  the  material  used  in  the  test 
h  <•  thickness  of  model 
n  =<  is  ochromatic  fringe  order 

Inasmuch  as  the  radius  equals  a  at  the  internal  boundary,  the  radial  stress,  (or  denoting  it  by 
the  principal  stress,  soy  <72)»  will  always  be  known:  c,  =  -p/  if  loaded  by  internal  pressure,  or 
Or  «:  0  if  unloaded.  The  maximum  (circumferential)  stress,  affiotO]),  will  be  completely  deter¬ 
mined  from  Eq.  1  as 

oa  -  2fn/h  +  Of  (at  the  star  point) 

providing  the  fringe  order,  n,  and  the  model  characteristics  are  known.  For  zero  internal  pressure, 
this  can  be  written  in  an  equivalent  form  as^ 


ae  _  Zf/h 
po  dpo/dn 


(2) 


In  the  problem  of  practical  interest,  that  is,  a  rocket  motor  subjected  to  internal  pressurization 
the  grain  is  subjected  to  lioth  an  internal  and  an  external  pressure  imposed  by  the  rigidity  of  the 
case.  As  will  bo  discussed  later,  these  continuity  and  displacement  conditions  at  the  case  are 
simulated  by  an  equivalent,  uniform  interface  pressure,  which  is  sufficiently  accurate  for  most 
engineering  purposes.  With  this  approximation,  the  important  parameter  in  determining  the  stress 
distribution  by  photoelustic  means  is  only  the  pressure  difference  across  the  web  (Ref.  1),  and  it 
may  bo  shown  ^  that  the  fringe  pattern  is  a  function  only  of  this  difference. 


CALIBRATION  TESTS 

For  the  calibration  tests,  it  was  necessary  to  determine  the  relationship  between  fringe  order 
and  stress.  Also,  since  the  specimen  was  to  bo  loaded  by  a  uniform  pressure  applied  to  the  outer 
circumference,  it  was  necessary  to  calibrate  the  loading  device.  The  device,  based  upon  tliat 
used  curlier  (Ref.  1),  applied  air  pressure  through  a  rubber  diaphragm  to  the  specimen.  Thus,  the 
relationship  between  air  pressure  indicated  in  a  gage  (pg)  and  the  pressure  applied  to  the  speci¬ 
men  (po)  was  required. 

The  fringe  constant  was  determined  by  loading  several  disks  of  various  diameters  with  two 
equal  diametral  loads  (Fig.  3).  The  load  required  to  produce  1,  2,  •••,  8  fringes  at  the  center  of 
each  disk  was  noted,  and  the  slope  of  the  straight  line  through  the  eight  data  points,  dL/dn,  was 
determined.  The  theoretical  solution  for  this  loading  is 


1  For  a  linear  system,  the  stress  Is  proportional  to  the  applied  pressure,  hence  og/po  =  Aag/dp^ 

=  (2//h)(rfn/dpo)  under  the  nssiintplion  that  at  zero  pressure,  the  grain  is  stress-free, 

A  uniform  pressure,  p,  applied  to  both  inner  and  outer  boundaries  of  any  multiply  connected,  two- 
dimensional  region  results  in  a  state  of  hydrostatic  stress:  i.e.,  Of  —p,  so  the  resultant  fringe  order, 

which  depends  on  0,-0$,  is  always  zero.  Thus,  this  state  may  be  superimposed  on  any  given  state  without 
altering  the  fringe  pattern.  A  problem  with  pt  internal  and  po  external  pressure  may  be  converted  to  a  problem 
with  Pi  —  Pi  =  0  internal  pressure  and  po  —Pi  external  pressure  without  changing  the  fringe  pattern.  Hence,  the 
fringe  pattern  depends  only  on  po  —  p<. 
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RECORDED  FORCE,  LB 

FIG.  3)  Typical  Slope  for  the  Determination  of  Fringe  Conatunt  for  CR-39  Plastic  Disk, 


dL  ^ 

dial -02)  '  8~  (3) 

where  L  is  the  load,  h  the  thickness,  D  the  diameter,  and  aj  ■~a2  is  the  difference  in  principal 
stress.  Thus,  the  value  ot2f/h  s  F  in  the  stress-optic  law  (Eq.  1)  ai  -02“  Fn  can  be  computed  for 

„  8idL/dn) 

Assuming  F  to  be  constant,  Eq.  3  becomes 

^  D  (4) 

dn  8 

The  average  value  obtained  from  the  tests  was  F  =  333  psi  per  fringe,  based  on  a  specimen  thick¬ 
ness  h  ®  0.250  inch.  Taking  the  value  of  E  to  be  333,  the  test  points  can  be  compared  to  the 
theoretical  value.  This  comparison  is  shown  in  Fig,  4,  and  reveals  only  a  small  scatter  of  data. 


dL/dn,  LB /FRINGE 
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FIC<  4<  Comparison  of  Tlieoretical  and  Experimental  Relationships  of  Load-to-Fringe 
Ratio  With  Diameter  of  Dlametrally  Loaded  CR-39  Plastic  Disics, 
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—  A  series  of  thick-^lled  Kdllw  cylindeiB  was  teatedjieja_ia^ie^ppgsaiire«loading  jig  to  deter- 
mine  the  relatlonsKip  beWeen^g  and  po*  ^  was  found  thatpg/po  ■  1?  hence,  no  further  distinction 
was  made  between  these  subscripts.  The  difference  in  principal  stress  in  a  thick-walled  hollow 
cylinder  at  the  inner  wall  under  external  pressure  loading,  po,  is  given  by 


ijlioi  —  O'^  2" 

dpo  1  -  Wb)^ 


- (5) 


Assuming  F  =  333,  dpo/dn  was  caiculated  from  the  test  data  and  compared  to  the  theoretical  value 
in  Fig.  5.  It  should  be  noted  that  the  data  fall  within  the  ±5%  band. 


FIG.  5.  Comparison  of  Theoreticai  and  Experimental  Relationships  of  Pressurc-to-Fringe 
Ratio  With  (t/b  for  an  Externally  Pressurized  Thick-Walled  Hollow  Cylinder. 
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STAR-CONFIGURATION  TESTS  ^ 

The  grain  configurations  tested  for  diis  report  were  all  of  the  type  shown  in  Fig.  1.  It  is 
clear  that  the  configuration  is  completely  specified  by  the  nondimenoional  parameters  w/b,  w/p, 
a’/a,  and  /V,  where  u>  ^  b  -  a,  and  that  the  stress  pattern  and  nondimensional  ratios  such  as 

we  independent  of  die  specific  magnitude  of  irFor  the  specimensTested^  6^^7480  inchesy 

andzoV^a  =  0rl37T^rhe  dependence  of  ffg/po  on  m/b,  w/prand  A'  is-presented  in  Figr6-T0r-An - 

empirical  equation  representing  the  data  appears  in  the  Appendix.  The  dependence  of  oe/po  on 
a'/o  was  investigated  for  two  configurations,  and  a$/pQ  was  found  to  be  nearly  independent  of  this 
parameter  for  values  of  a'/a  up  to  0.6.  These  results  are  presented  in  Fig.  11. 


FIG.  64  Effect  of  Grain-Dimension  Parameters  on  the  Stress  at  the  Star  Point  of  an  Ex¬ 
ternally  Pressurized  Three-Poinl-Star  Grain, 


The  specimens  to  be  tested  were  machined  from  CR-39  plastic  on  a  milling  machine.  The 
procedure  was  such  that  the  dimensions  of  the  star  configuration  were  held  within  +0,005  inch, 
and  the  angular  position  of  tlie  star  points  was  held  within  ±5  minutes  of  arc.  These  tolerances 
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FIG.  10.  Effect  of  Grain-Dimension  Parameters  on  the  Stress  at  the  Star  Point  of  an  Ex¬ 
ternally  Pressurized  Eight-Pcint-Star  Crain. 
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0  0.2  0.4  0.6  0.8  1.0 


o' fa 

(n)  Thrae-polnt-slar  grain. 


(b)  Four-point-star  grain. 

FIG.  11.  Effect  of  Inverse  .Star  Point  on  the  Stress  at  the  Star  Point!  o/h  =  0,5j  tc/p  =  10.79. 


are  more  than  enough  to  ensure  that  any  errors  in  the  final  result  arising  from  inaccurate  geometry 
are  completely  negligible  compared  to  other  errors  inherent  in  the  method. 

CR-39  exhibits  a  time-edge  effect  that  tends  to  obscure  and  distort  the  fringes  near  the  edge 
of  a  specimen  when  it  is  allowed  to  stand  after  machining.  A  lapse  of  over  an  hour  is  enough 
time  to  hamper  severely  a  photoelastic  determination  of  edge  stresses;  therefore,  each  configura¬ 
tion  was  tested  immediately  after  machining.  By  widening  and  extending  the  star  point  at  each 
step,  it  is  possible  to  obtain  several  specimens  cut  to  different  configurations  from  one  blank, 
while  producing  a  fresh  edge  to  minimize  time-edge  effect  at  each  step.  In  this  way,  from  three 
to  seven  test  points  may  be  obtained  from  each  blank. 

As  a  control,  each  blank  was  first  machined  to  a  thick-walled  cylinder  configuration  with 
an  inside  diameter  of  0.375  inch  (corresponding  to  a/h  »  0.137)  and  pressure-tested.  After  com¬ 
pletion  of  the  parametric  changes  of  star-configuration  tests  on  a  given  blank,  the  remainder 
was  again  cut  into  a  thick-walled  cylinder  and  tested.  These  two  tests  of  a  configuration  whose 
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stresses  could  be  analvzed  theoretically  pave  two  check  points  on  each  blank^IiradditioD;  a 
small  disk  was  cut  from  the  remains  of  each  blank  and  tested  by  applyiJig  a  diametral  load>s 
previously  described.  Since  the  two  loading  rigs  were  different,  this  procedure  provided  an  in¬ 
dependent  check  of  the  fringe  constant  for  each  blank. 

_ llEFFECT-OF  THE  INVERSE  STAR  POINT 

Another  aeries  of  tests  was  run  to  determine  the  effect  of  the  shape  of  the  inverse  star  point 
(point  A  of  Fig.  1)  on  tlie  stress  at  the  star  point.  The  purpose  of  determining  this  effect  was  two¬ 
fold!  (1)  so  the  designer  would  know  whether  the  results  presented  herein  were  applicabJe  to  the 
grain  in  question  and  (2)  to  give  him  an  indication  of  the  effect  on  the  stress  at  the  star  point 
caused  by  minor  variations  of  the  internal  geometry  that  might  be  imposed  for  ballistic  reasons. 

Two  types  of  tests  were  conducted.  The  first  and  most  significant  was  to  vary  a' while  the 
remainder  of  the  geometry  was  held  constant.  The  initial  point  on  each  curve  is  at  a  value  of 
a'/a  thnt  corresponds  to  the  pointed  inverse  star  point.  This  condition  has  the  relationship 

a"  p  rr 
—  -  —  CSC  — 

«  rt  N 

(Smaller  values  of  a/a  do  not  correspond  to  u  possible  star  configuration.) 

Tho  tests  wore  made  liy  iiiachiniiig  a  blank  to  tlie  pointed  inverse  star  configuration,  testing  it 
while  it  was  being  subjected  to  cxtenuil  pressure,  and  then  rcm.icliiiiing  it  to  increase  a'/a,  'I'ho 
apucimun  vvas  tested  and  the  process  rc|ieate(|,  liecausc  the  configuration  remained  tho  same  in 
the  region  near  the  star  point  where  the  inaxinuim  stress  was  being  detcrinined,  a  frosli  edge  was 
not  produced  there  before  each  test;  lienee,  the  ti<ue«edge  effect  must  be  accounted  for,  Tho  solid 
curves  of  Pig.  II  are  based  on  the  iiiicon'cctud  tost  data;  the  dashed  curves  represent  data  that 
liavo  been  corrected  for  tiiuo-odgc  effect  by  using  some  recently  obtained  (tontntive)  correction 
factors,  At  this  time,  tlioro  are  insufficient  data  to  decide  between  tliese  curves,  IJow'ever,  the 
true  effect  is  undoubtedly  bracketed  by  the  corroded  and  iincorrocted  curves. 

The  second  type  of  test  vvas  run  to  check  results  presented  by  Wilson  (liof,  3  und  4)  for  an 
analytic  solation  that  will  be  discussed  later  in  lliis  report,  and  to  check  the  effect  of  changing 
the  entire  cliuracter  of  the  internal  boundary.  In  these  tests,  the  internal  boundary  was  made  to 
correspond  to  Wilson’s  geometry  os  closely  ns  conveniently  possible.  This  resulted  in  a  square 
hole  rather  than  the  normally  circular  one.  1'he  results  of  these  tests  were  not  significantly  dif¬ 
ferent  from  those  of  an  inscribed  circubu  hole. 


PROCESSING  OF  TEST  DATA 

The  pressure-loading  tests  were  made  by  noting  the  pressure  at  which  each  fringe  appeared  at 
the  point  where  stress  was  to  be  determined.  In  the  thich-walled  hollow  cylinder,  the  fringes  were 
rather  broad  rings,  which  made  an  accurate  determination  of  the  location  of  the  fringe  difficult.  In 
most  cases,  only  four  fringes  could  be  observed  because  of  a  400-psi  equipment  limitation  on  the 
pressure.  I  he  combination  of  indistinct  fringes  and  four  or  five  data  points  makes  the  deter¬ 
mination  oidpo/dn  (and  thus  da^/dpo)  tlifficult.  Typical  data  points  can  be  observed  in  Fig.  12. 

In  contrast,  the  star  configuration  exhibits  a  very  sharp  fringe  pattern;  eight  to  ten  fringes  can  be 
observed  without  exceeding  pg  =  400  psi.  Also,  these  eight  to  ten  data  points  were  obtained  for 
at  least  three  star  points  on  every  configuration,  giving  three  sets  of  test  points  from  which  the 
slope  dpo/dn  could  be  determined.  Typical  data  are  shown  in  Fig.  13. 
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FIG.  12.  FrinBo  Order  Vothus  Pressure  for  un  Kxtcrnnlly  Pressurized  Thick-Wnllod 
Hollow  Cylinder  of  CH-39  Pliistici  a/h  =  0.368;  h  =  0.253. 
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FIG.  13.  Fringe-Order-Versus-Pressure  Curves  for  Three  Star  Points  (A.B.C)  of  an  Ex 
temally  Pressurized  Six-Point*Star  Grain. 
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The  dataJrojrLeactlesL  were  plotte^d  in  the  po  versus^  B-pJenef  anitlre  best  straight  line - 

(visuallyl  was  plotted  througbrthe  datarrThe  slope  of  tbis  line  was  read  fronrthe^lotT^A^MlMi 
lation  method  was  also  used  that  determined  the  slope  of  the^traight  line  thrdugh^he^ata  points 
by  minimizing  the  sum  of  the  squares  of  the  distance  of  the  data  points  from  the  line  (Ref.  5). 

The  slope  determined  in  this  manner  was  compared  to  the  visually  determined  slope.  It  soon 

became  apparent  that  the  leaat-jquare^calculated  slopes  yielded  data  of  slightly  greater  jilternaL _ 

consistency.  The  scheme  also  eliminated  the^  human  element  and  provided  a  standard  way_of_de- 
fihiTig  the  “best  slope.*  Ultimately,  the  plotted  data  became  simply  a  chech  against  wild  (mis¬ 
read)  data  points  or  calculation  arrors.  The  calculated  values  of  slope  were  used  for  all  subse¬ 
quent  plotting  and  calculation  and  are  denoted  by  (dpo/dn)calcd. 

In  view  of  the  more  accurately  determined  slopes  dpo/dn  for  the  star  configuration,  it  is  clear 
that  the  curves  determined  from  these  data  must  compare  to  the  unknown  theoretical  values  at 
least  as  well  as  the  poorer  thick-walled-cylinder  calibration  data  compare  to  the  known  theoreti¬ 
cal  curves,  as  shown  in  Fig.  .5. 

A  total  of  .34  data  points  on  the  doo/dpo  versus  tv/p  and  w/b  curves  was  obtained  for  the  six- 
point-star  configuration,  which  was  chosen  us  the  control  on  which  the  most  extensive  work  would 
be  done.  In  addition,  a  limit  point  was  known  on  each  curve  of  constant  w/b  when  plotted  versus 
w/p.  Consider  the  situation  when  w  and  h  are  held  constant  and  w/p-*(^.  It  is  clear  that  this  process 
implies  p->oo.  However,  a  moment’s  reflection  on  tlio  geometry  indicates  that  the  star  configuration 
ceases  to  exist  for  p  >  a  and  that  the  case  p  r;  a  reduces  to  u  thick-walled  cylinder.  (Jd/po  for  this 
cylinder  cun  he  ciilciilated  from  hlcj.  .5,  and  represents  the  limit  point  (shown  ns  crosshntclied  ureas 
in  Fig.  (>—10)  on  the  curves  as  u'/p-»0.  With  these  limit  points  pins  four  experimental  points  on 
each  curve,  the  shapes  of  the  curves  of  constant  tv/b  in  the  o$/pQ  versus  iv/p  plane  were  cpiito 
wel I  defined. 

For  the  three-,  four-,  five-,  and  eight-point  stur  configurations,  a  total  of  .30  tests  was  run. 

The  limit  point,  tlin  two  dat*i  points  on  ouch  of  these  curves,  nnd  the  known  shape  of  the  curve 
from  the  six-point  star  data  were  siifficicnt  to  define  the  curves  for  these  configurations.  Cross¬ 
plots  of  various  purainotors  versus  A'  were  also  constructed  from  the  initial  curves  and  used  to 
smooth  the  Og/po  versus  w/p  curves.  An  example  of  a  crossplot  is  given  in  Fig,  14. 

ERROR  IN  EXPERIMENTAL  RESULTS 

The  final  results  are  given  in  the  funne^/po  versus  w/p  along  curves  of  constant  w/h,  'I’he 
location  of  a  given  experimental  data  point  in  the  og/pQ  versus  w/p  plane  is  determined  by  these 
three  values. 

As  bus  been  previously  pointed  out,  mucliining  errors  resulting  in  inaccurate  geometry  are 
entirely  negligible;  therefore,  the  point  may  be  assumed  to  be  perfectly  located  with  respect  to 
the  two  parameters  that  depend  only  on  geometry  (w/p,  w/b).  Thus,  the  error  of  location  of  an 
experimental  point  reduces  to  an  evaluation  of  the  error  in  deleritiining  ag/po.  In  Kq.  2,  the  rela¬ 
tive  error  in  ag/pQ  is  the  sum  of  the  relative  errors  in  2f/h  =  F(0.250)/A  and  in  dpa/dn.  From 
Fig.  4  it  may  be  seen  that  the  error  in  F  is  within  ±0,5%.  3'he  thickness  of  each  specimen,  h, 
was  measured  with  a  micrometer.  The  thickness  was  nominally  0.2,30  inch,  and  readings  were 
taken  to  the  nearest  0.001  inch.  However,  the  thickness  tended  to  vary  over  a  given  specimen, 
so  the  average  value  was  used  in  calculations.  Specimens  with  a  variation  of  more  than  +0.003 
inch  were  rejected.  Using  the  average  value  of  h  under  these  circumstances  certainly  could  not 
introduce  an  error  of  over  0.003/0.250  x  100%  =  1,2%,  However,  the  expected  error  would  be  con¬ 
siderably  less  than  this,  and  +0,.5%  represents  a  conservative  estimate.  Thus,  the  total  contribu¬ 
tion  to  error  from  f/h  is  +1%. 
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readlt‘'!'‘“rrTT  ^p.Un  is  mon.  difficult  to  estimate.  The  pressure  was 

read  from  a  calibrated  gage,  so  the  error  from  this  source  is  neRliKible  Eiaht  to  ton  dof 

W...  u«,d  to  d„o™|„  „ob  .lmight.|i.,  curv.,  .od  thro.  i«d.pojdo„t\.tjo.  ^te  „ben /Z 
..oh  Th,  ,t  the  daa  i.  otoelle;,  (Fig.  U) 

A  maximum  of  ±3%  relative  error  in  the  (averaged  over  three-star-noint)  ,/n  /,/  d  . 

".ll.d.0,li„d»  d...  (L4.~  Fig  12  ““'“•'“I'  o'  >i"  thioh- 

boslIlZiZ''  "’’"“'’"‘f'  ■‘o  ""•''0.  .hoold  b.  ,lthi„  ±4*  otth.  cortoM  The 

b..t  .i„o.lh.d  curve.  „i,|,  lb,  ,  i„die„.d,  i„  Fig  6-10  iT  J!’ 

ere  .....»d  b.  ooutinuou.,  »„.,.„„ie,  ..d  ...oZbl,  ...,.h  „i,hZ^cUo  .iu!‘ 
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- then^he  final  curves  probably  fall  well  within  ±4%  of  the  over  Wst  orthe^^^^ 

engineering  interest, - - - 


COWPARISON  WITH  A  THEORETICAL  SOLUTION 


Jil8oiL(H^.  J^and  4)  determined  the  stress  in  a  four-point-star  grain  analytically  bywaw  of 

a  complex-variable  approach.  In  this  analysis,  he  uses  a  four-term  transformation  that  transforms  a 
circle  into  a  shape  matching  the  shape  of  the  grains  tested  herein  in  two  parametets!  namely,  po¬ 
sitions  of  star  points  and  radius  of  curvature  at  the  star  point.  The  over-all  shape,  however,  is 
somewhat  different,  particularly  in  regions  away  from  the  star  point.  Fortunately,  the  effect  of 
these  regions  on  the  maximum  stress  at  the  star  point  is  negligible  (Fig.  11). 

The  results  of  the  experimental  test  for  a  configuration  similar  to  that  of  Ref.  4  are  shown  in 
t  ig.  15,  which  gives  a  direct  comparison  of  the  two  solutions  over  the  common  range.  As  may  be 
seen,  the  agreement  is  excellent;  the  volues  of  the  test  data  ore  at  most  2%  higher  than  the 
theoretical.  ” 


EXTENSION  OF  TEST  DATA 
STRESS-CONCENTRATION  FACTORS 

Although  the  results  presented  thus  far  are  for  a  grain  subjected  only  to  external  pressure, 
the  results  may  bo  extended  to  several  other  cases  of  practical  interest. 

Fr*#  Grain  With  Internal  Pretsure.  For  a  groin  that  is  unloaded  at  the  outer  boundary  but  sub- 
jeeted  to  an  internal  pressure,  p{,  the  use  of  Kq.  1  and  2  putting  =  og  and  aj  =  -p,'  gives 


HG.  15.  Comparison  of  Theoretical  and  Experimental  Stress-Concentration  Factors  for 
an  Externally  Pressurized  Four-Poinl-Slar  Grain;  a/b  =  0.5. 
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a$ - 2f/h 


-1 


(6) 


P(  dpi/dn 

or,  because  (2f/h)/{dpi/dn)  is  formally  equivalent  to  the  value  a$/pQ  in  Fig.  6—10, 

Vfl  as  ^ 

_  Pt  Po 

Hence,  the  stress  at  the  star  point  eaused  by  an  internal  pressure,  pft  may  be  determined  immedi¬ 
ately  by  reading  the  value  of  as/p^  off  the  graphs  at  the  appropriate  values  of  the  geometric  param¬ 
eters  and  using  Eq.  6, 

CatS'Bended  Grain  With  Internal  Pressure.  For  a  case-bonded  configuration,  it  is  convenient, 
although  not  necessary,  to  introduce  the  concentration  factor  defined  in  Ref.  2  as 


-  ffj 

= ::: - 

CTj  -  Oj 


(7) 


where  the  barred  quantities  refer  to  the  stress  in  a  corresponding  thick-walled  cylinder.  For  the 
case  in  question. 


-  _  _  2(po-p,.) 

Oi  -  Oj  =0,-00^  - — •— 


(8) 


1  - 

rile  factor  2/[l  -{a/b)^\  has  boon  plotted  versus  w/b  in  Fig.  16.  Tlierefore,  the  concentration 
factor  may  bo  determined  immediately  from  the  expression 

y  _  fox 

'  2/tl-(a/A)^] 

The  value  of  as/po  •«  obtained  from  Fig.  6-10.  It  is  of  interest  to  note  timt  the  expression 

(10) 

holds  for  all  the  expressions  that  follow.  Hence,  ns  indicated  above,  the  introduction  of  the  stress- 
concentration  factor  is  not  necessary,  but  was  done  so  that  results  quoted  herein  might  be  com¬ 
pared  with  earlier  data. 

In  actual  practice,  the  grain  will  he  restrained  by  a  relatively  rigid  case.  As  stated  previously, 
it  has  been  found  satisfactory  for  engineering  accuracy  to  replace  the  effect  of  the  case  enclosing 
an  internally  pressurized  star  grain  by  p',  a  uniformly  distributed  interface  pressure.^  If  the  thick¬ 
ness  of  the  case  is  assumed  small,  the  interface  pressure  is  given  by 

2p( 

p'  - ;Vv2-77rT2 — ?r- 


1  + + 


2(1  -  y)pi 


‘  .1 


e.K 


1  + 


O  i-b 

(1  -2(/)  -  +  -  -1  - ^ - 

\o/  J  l\o/  J  1  + 1/ 


(plane  strain) 


(12) 


3  The  actual  interface  pressure  caused  by  internal  pressurization  varies  because  of  the  irregular  shape 
of  the  star  grain.  This  variation  is  well  approximated  by  a  sinusoidally  varying  pressure  applied  at  the  outer 
boundary.  However,  the  authors  in  unpublished  experiments  have  found  that  the  effect  of  this  variation  on  the 
stress  at  the  star  ooint  is  small* 
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h-/«,we0  fraction 

riG.  I6«  Vnrifltlon  of  Normalizing  Factor  With  Wol)  Fraction. 


where  /t^  is  the  thickness  of  the  case  and  the  subscript  c  on  other  quantities  refers  to  the  case. 
From  Eq,  7, 

<7,  -  -  as)  (13) 


The  normalizing  stress  is  given  by  Eq.  8  and  11  or  12,  whichever  is  applic,  ' 
at  the  inner  surface  is  -pi.  Inserting  these  values  in  Eq.  1.3,  and  for  simplicity  ass 
compressible  grain  (iz  -  1/2)  and  a  state  of  plane  strain,  the  expression  for  the  tangential  stress 
becomes 


2Ki 


1 


(14) 


or,  using  Eq.  10, 


Pi 


1 


(15) 


Hence,  the  stress  under  such  a  loading  condition  is  determined  from  the  expressions  above,  and 
the  value  of  ff^/po  is  obtained  from  Fig.  6—10. 
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STRAIN^CONCENtRAflON  FACTORS^ _ 

Strain-concentration  factors  may  also  be  calculated  from  the  information  contained  in  this 
report.  The  stress— strain  relationships  in  cylindrical  coordinates  for  a  state  of  plane  strains  are 
given  by 

fr  =  (16) 

=  a(ag  +  /3a,)  (17) 

y,B  ^  (1/G)r,e  (18) 

where 


a  = 


1 


V 

iTTi 


(19) 


A  strain-concentration  factor  may  be  defined  in  a  manner  similar  to  the  stress-concentration 
factor,  i.e., 


<1  -  f2 


(20) 


whore  and  <2  arc  the  principal  strains,  and  the  barred  quantities  denote  the  strains  in  a  thick- 
wnllod  hollow  cylinder.  In  an  isotropic  elastic  material,  the  directions  of  principal  stress  and 
strain  coincide;  hence,  along  the  lino  of  symmetry  where  there  is  no  shearing  strain,  Kq.  20  may 
be  written  as 


,,  (r  -  ie 

K  ~ - :r- 

it  -  *6 


(21) 


Using  fclq.  16  and  17,  the  difference  in  principal  strain  may  ho  written  us 

(t'~  i3  -  “  l)(oi>  “  a,)  (22) 

Hence,  ns  would  bo  expected,  the  difference  in  principal  strains  is  n  function  only  of  the  differ¬ 
ence  in  principal  stress,  which  in  turn  depends  only  on  pressure  difference.  This  fact  allows  any 
stotements  or  proofs  concerning  the  relationship  of  the  fringe  pattern  resulting  from  internal  or  ex¬ 
ternal  pressure  in  the  stress  problem  to  be  applied  directly  to  the  strain  problem. 

Now,  for  the  specific  problem  under  consideration:  namely,  internal  pressurization  of  n  case- 
bonded  grain, 

a,  =  -Pi 

2Ki 


ag  = 


1  -  (o/i) 


;(1-  P)ipi)-P( 


(23) 


a(/9 


p.t 

Pi 

Fur  the  normalizing  cylinder,  all  quantities  are  the  same  with  the  exception  of  the  concentration 
factor,  which  is  unity.  Combining  the  expressions  of  Elq.  23  in  Fq.  22  gives 

1+1/  2K( 


=  • 


l-(a/fc)2 


iP-l)pi 


(24) 
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hence 


Kt  =  =  Ki  (25) 

(,  -  (e 

That  IS)  the  nunierical  value  of  the  strain-concentrution  factor  as  defined  ^bove^or  anzinternally 
pressurized  grain  is  identical  to  the  stregs-cpncentration  factotiar  e  free  grain  subjected  to  ex¬ 
ternal  pressure. 

It  is  important  for  the  designer  to  note  that  the  radiaj  component  of  the  strain  does  not  vanish 
at  the  star  point.  In  the  case  of  the  stress-concentration  factor,  the  radial  stress  was  either  zero 
or  the  negative  of  the  internally  applied  pressure;  in  an^  event,  the  radial  stress  w^as  specified 
independently  of  the  tangential  stress.  Elxatnination  of  Eq.  16  reveals  that  this  is  not  true  for  the 
strain-concentration  factor.  There  is  a  contribution  to  the  radial  strain  from  both  the  tangential 
and  radial  stresses.  It  is  possible  that  the  tangentiol  strain  or  the  ratio  of  the  tangential  strains 
to  the  normalizing  tangential  strains  might  become  quite  large  compared  to  the  strain-concentration 
factor  defined  above.  Hence,  the  designer  must  decide  which  parameter  is  important  in  the  failure 
of  a  grain  and  compute  the  concentration  factor  from  the  data  presented  accordingly.  The  funda¬ 
mental  importance  of  establishing  an  appropriate  failure  criterion  has  been  strongly  emphasized 
(Hof.  2). 

'Iho  extensions  pi'esontod  above  do  not  exhaust  the  list  of  possible  ones  that  can  be  made 
from  the  data  given  heroin.  Indeed,  extensions  to  certain  thermal-stress  situations  have  been 
pointed  out  (Hof.  1  and  6).  Finally,  those  data  do  not  cover  all  configurations  of  ballistic  interest, 
but  it  is  believed  that  tliis  pnramotric  study  of  one  set  of  typical  geometries  will  help  the  designer 
at  least  to  estimate  the  factors  for  his  specific  problem. 

CONCLUSIONS 

Photoolastic  tests  wore  conducted  to  determine  the  maximum  tangential  stress  for  a  family  of 
star-grain  configurations,  (lurves  are  presented  that  give  the  value  of  this  parameter  for  web  frac¬ 
tions  from  20  to  60%,  for  w/p  from  approximately  2  to  30  (where  w  is  the  web  thickness  of  the 
groin,  and  p  is  the  radius  of  the  star  point),  and  for  N  (number  of  star  points)  from  3  to  8.  The 
curves  are  expected  to  be  accurate  to  ±4%  over  the  range  tested. 

The  effect  of  the  inverse  slur  point  was  investigated  and  found  to  bo  negligible  over  a  rongo 
of  shapes  and  values  of  the  parameter  a'/altom  approximately  zero  to  0.6,  The  experimental  re¬ 
sults  are  compared  to  one  known  theoretical  calculation,  and  the  agreement  is  found  to  be  excel¬ 
lent. 
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Appendix 

DERIVATION  OF  EMPIRICAL  EQUATION 

The  results  of  the  tests  described  in  the  body  of  this  report  are  presented  in  Fig.  6—10.  It  is 
the  purpose  of  this  appendix  to  show  the  derivation  of  an  empirical  equation  that  represents  these 
experimental  data.  The  equation  derived  gives  the  ratio  of  tangential  stress  to  external  pressure 
as  on  explicit  function  of  the  three  controlling  parameters:  web  fraction,  radius  of  the  star  point, 
and  number  of  star  points.  The  ranges  of  these  parameters  are  as  follows:  web  fractions,  from  20 
to  60%;  w/p,  from  2  to  30  (where  w  is  the  web  thickness  of  the  grain  and  p  is  the  radius  of  the 
star  point);  and  A'  (number  of  star  points),  from  3  to  8.  The  accuracy  of  the  equation  was  demon¬ 
strated,  and  good  agreement  with  the  experimental  data  was  achieved. 

The  main  reason  for  obtaining  an  explicit  algebraic  relationship  for  the  stress  ns  n  function 
of  the  various  paramoters  is  tiint  such  a  relationship  may  he  used  directly  in  ballistic  calculations 
Previously  it  was  necessary  to  arrive  at  a  configuration  for  n  specific  ballistic  requirement  by  an 
iterative  process.  In  such  a  process,  a  geometry  is  first  selected  without  regard  to  the  stress 
level  involved.  This  configuration  is  then  tested,  and  the  results  obtained  arc  used  to  judge  the 
relative  merit  of  the  design  from  the  stress  standpoint,  ilowcvcr,  with  the  equation  presented  and 
by  the  use  of  n  high-speed  computer,  it  is  possible  to  introduce  stress  considorntions  directly  into 
the  selection  process.  In  principle,  for  n  given  ballistic  mission,  the  minimum  stress  configurotion 
of  n  family  of  grains  could  be  doteririinod  purely  by  ninchine  computation. 

METHOD  OF  CURVE  FITTING 

The  nature  of  the  oxperimoiitnl  curves  presented  in  Fig.  6—10  and  crussplots  such  ns  Fig,  H 
suggest  that  the  data  may  be  represented  by  a  function  of  the  form 

op/po  =  A'*'d»(M'/i,  lo/p)  (26) 

If  w/b  and  w/p  are  hold  constant  (i.c.,  the  function  ^),  the  value  of  4:  may  be  obtained  from  the 
experimental  curves,  Uy  the  application  of  this  procedure,  the  assumption  that  /c  is  independent  of 
w/b  and  w/p  was  verified,  except  in  a  limited  region  near  the  limit  point 

Ml  w/b 
p  1  -  (w/b) 

The  value  of  /e  thus  determined  is  -0,34. 

Knowing  the  dependence  of  the  curves  on  A',  the  functional  nature  of  was  determined  from 
the  curves  for  a  specific  value  of  A'.  The  curves  for  A'  -  6  were  chosen  for  this  purpose.  The 
shape  of  these  curves  suggested  choosing 

<t>(w/b,w/p)  =  a.{w/b).(w/p)^'^'’’K  y(w/b)  (27) 

or 

log(4i-  y)  =  /3  lo^(w/p)  +  log  (X  (28) 

For  fixed  w/b('20%,  30%,  •  •  •  ,  60%),  lOq.  28  gives  the  relationship  between  the  experimentally 
known  ft>(w/p)  and  the  constants  U  and  /3,  The  right-hand  side  is  the  equation  of  a  straight  line 


20 


_ ii-tlLeJog  log  plane,  vyith  slope  /3  aiid-intorcept  ai-By-ohoosing  several  values  of  yrit-wasTpos^ 

sible  to  raake-the  experimentally  determined  data  (left»hand  side)  also  plonsa^traighrliSe^in  tte 
log  log  plane,  thus  showing  that  the  assumed  form  for  the  function  <&  is  an  appropriate  one.  The  y 
that  makes  the  experimental  data  plot  as  a  straight  line  and  the  CC  and  /3  determined  from  this 
straight  lino  were  then  known  for  each  value  of  w/b,  Ot^jSj^andj/ were  plotted  as  a  function  of 
w/bf  and  appropriate  functions  were  fitted  to  the  cwws^o^ defined.  In  this  manner,  the  following 
functions  were  determined 

a  =  0.52  Wh) 0.80 

/3  =  1.90  1.0  (29) 

y  *  1.7(1  -  15.8  exp  [-18,4(m/6)]i 
The  empirical  function  thus  obtained  is 

(Te/po  =  (6/A')“-3*([0.52(«./i)-»-77  _  0.80](,i;/p)  Cl.90(«./6)0-‘75-i,o] 

+  1.711  -  15.8  exp  N8.4(m/W]|)  (30) 

This  function  holds  over  the  region  of  the  experimental  data,  hut  fails  near  the  limit  point 

If  w/  b 

p  1  -  {w/b) 

Near  this  point,  tlio  value  of  og/po  nppmnehos  the  value  for  n  tluelc-wallecl  cylinder,  and,  thus, 
cp/po  is  independent  of  A',  To  extend  the  validity  of  Eq,  30  to  this  limit  region,  tlie  term  describ¬ 
ing  tlio  A'  (Icpondcnco  was  modified  in  llic  following  manner!  n  function  was  devised  that  tended 
to  (6/A')“‘3^  for 


W  w/b 

-  » - 

p  i  -  iw/h) 

and  tended  to  1  ns  w/p  tended  to 

w/b 

T  -  iw/b) 

A  satisfactory  function,  which  provides  n  smooth  transition  for  the  N  dependence  of  this  nature. 


If  w  w/b  b 

-114 1 - ; -  - 

Ip  1  -(w/b)Jw. 


Henco,  in  any  problem  in  which  the  value  of  ap/po  near  the  limit  point  is  required,  the  following 
expressions  should  he  used  instead  of  Eq.  30. 


»  1  + 


f/6r  1/  f  f*"  »<v 

l(/v)  “fiTi 

/w\-1.77 


w/b  b 
■  {w/b)jw 


{[0.52(j)-‘-"-0.8oJ(*) 


+  1.7  1  -  15.8  exp 


('»■%- )ll 


It  should  be  noted  that  for  normal  usage,  the  range  covered  by  Eq.  30  will  be  sufficient. 
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ACCUI?ACY^F^RKUUrS 

Values  of  O's/po  given  by  Cq.  32  in  the  region  of  the  experimental  data  were  calculated. 
Curves  were  drawn  through  these  points  and  superimposed  on  the  experimental  curves.J^ryes 
for-the  three-pwit-^gtar  grain  are  showndnrFig.  17.  It  will-be  obsMvedTthat  the  agreement  is  ex- 
cellent.  with  discrepancies  of  less  than  5%  over  oLOst  of  the  Jange  of  interestr^he  maximum 
deviation  occurs  on  tlie  threespoint-star  curve  for  w/b  of  40%.  This  single: polfit  deviatej„ftQin 
the  experimental  one  by  11%;  It  should  be  recalled  that  the  accuracy  of  this  curve  was  quoted  as 
+  4%!  hence,  agreement  might  range  between  7  and  15%.  While  the  error  of  this  single  curve  at 
first  observation  appears  large,  it  must  be  viewed  in  thejight  of  the  excellentagreement  of  the 
empirical  equation  with  the  remainder  of  the  data  and  the  fact  that  this  curve  was  determined  by 
only  two  data  points.  This  discrepancy  could  easily  be  corrected  by  adding  yet  another  factor  to 
Eq.  32.  However,  because  of  the  over-all  agreement  with  the  remainder  of  the  data,  the  additional 
complication  Involved  was  not  considered  justified. 


FIG.  n.  Comparison  of  Experimental  and  Analytical  Curves  for  the 
Three-Point-Star  Grain. 
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